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ABSTRACT 

We estimate the rate of near-field microlensing events expected from all-sky surveys and investigate the 
properties of these events. Under the assumption that all lenses are composed of stars, our estimation of the 
event rate ranges from Ftot ^ 0.2 yr"' for a survey with a magnitude limit of = 12 to Ftot ^ 20 yr"' for a 
survey with Vnm =18. We find that the average distances to source stars and lenses vary considerably depending 
on the magnitude limit, while the dependencies of the average event time scale and lens-source transverse speed 
are weak and nearly negligible, respectively. We also find that the the average lens-source proper motion of 
events expected even from a survey with Vwm =18 would be (/i) > 40 mas yr~', implying that the source and 
lens of a significant fraction of near-field events could be resolved from high-resolution follow-up observations. 
From the investigation of the variation of the event characteristics depending on the position of the sky, we find 
that the average distances to source stars and lenses become shorter, the lens-source transverse speed increases, 
and the time scale becomes shorter as the the galactic latitude of the field increases. Due to the concentration 
of events near the galactic plane, we find that > 50% of events would be detected in the field with b < 20°. 
Subject headings: gravitational lensing 



1. INTRODUCTION 

The concept of one star gravitationally amplifying the 
light fro m another background star was first considered by 
lEinsteinl ([T936i) . although he concluded that the chance to ob- 
serve this phenomenon (gravitational microlensing) would be 
very low. With the advance in technology, however, detec- 
tions of microlensing events became feasible from system- 
atic se arches. The first detections were reported in the early 
1990s dAlcock et al.|[T993l:lAubourg et al.|[l99llUdalski et all 
Il993h . To date, lensing events are routinely detected with 
a rate of > 500 events per season and the total number of 
detections now exceeds 3000. To maximize the number 
of detections, these searches have been and are being con- 
ducted toward very dense star fields of nearby galaxies such 
as the Magellan ic Clouds jX lcock et al. 2000; Afonso et al] 
^03") and M31 (Ansai-iet an999; Calchi Novati et al. 2005^ 
Uglesich et al. 2004; de Jong et al. 2004; Riffeser et al. 2003), 
and galactic center dAlcock et al. 200 U IHamadache et all 
l2006tlSumi e"tani2006HBo'ndet aL,.200C 

Recently, a detection of a le nsing event o ccurre d on 
a nearby s ta r was reported by iGaudi et akl (l2007h and 
iFukuiet al.l ( IIoOTI) . The star (GSC 3656-1328) is lo- 
cated about 1 kpc from the sun in the disk of the Milky 
Way. This detection along with the advent of transient sur- 
veys capable of cove ring a very wide field with high ca- 
dence such as AS AS ('Szczvgiel& Fabrvckv"2007), ROTSE 
dAkerlof et al. 2003), TAROT and ARAGO (Boer 2001), and 
Pan-Starrs (Hodapp et al. 2004) have drawn attention of many 
researchers in microlensing community on the feasibility of 
near-field microlensing surveys. If near-field events could be 
detected from such surveys, they would be able to provide 
precious information about the matter distribution around the 
sun including dark objects. In addition, these surveys might 
enable detections of close planets with the microlensing tech- 
nique th at has demonstrated its capability in detecting distant 
planets ("Bond et al."2004l; IUdalski et alj|2005l; iBeaulieu et alj 
[2006; Gould et al. 2006). 

The probability of lensing occurring on nearby stars has 



been estimated bv'Collev & G ottl(ll995h an d lNemiroffl(ll998l) . 
The estimation of Collev & Gotll ( 1995h focused on only 
events that could be observable with naked eyes, i.e. with 
a magnitude limit Vum ^ 6. Nemiroff (1998) estimated the 
lensing probability expected from surveys with various mag- 
nitude limits. However, he estimated the probability in terms 
of the average number of stars undergoing lensing magnifica- 
tion at a moment, not in terms of events per year, by using the 
combination of the star-count information and optical depth 
to lensing. The optical depth to lensing represents the proba- 
bility that any given star is microlensed with a magnification 
A > 1 .34 at any given time. This optical depth is dependent on 
the mass distribution but it is independent of the mass of the 
individual lensing objects. Therefore, the probability based on 
the optical depth does not provide detailed information about 
the physical parameters of lens systems such as the lens mass, 
locations of the lens and source, and their relative transverse 
speed, and the observables of events such as the event time 
scale and source brightness. For this information, additional 
modelling of the mass function and velocity distribution of 
galactic matter is required. 

I n this paper , we ex tend the works of ICoUey & GottI (Il995h 
and iNemirofj (Il998l) to estimate the rate of lensing events 
in terms of event per year and to obtain the distributions of 
the physical parameters and observables of events expected 
from all-sky lensing surveys with various magnitude limits. 
For this, we conduct Bayesian simulation of near-field lens- 
ing events. 

The paper is organized as follows. In § 2, describe the de- 
tails of the simulation. In § 3, we present the resulting event 
rate and distributions of the physical parameters of lens sys- 
tems and observables of lensing events. We also investigate 
the variation of the event characteristics and lensing proba- 
bility depending on the position of the sky. We analyze the 
tendencies found in the distribution and explain the reasons 
for these tendencies. We summarize the result and conclude 
in §4. 

2. LENSING SIMULATION 



2 



NEAR-FIELD MICROLENSING 



To investigate the rate and properties of near-field lens- 
ing events expected to be detected from wide-field surveys, 
we conduct Bayesian simulation of these events. The basic 
scheme of the simulation is as follows. 

1 . We first produce source stars on the sky that can be seen 
from a survey with a given magnitude limit. We assign 
the locations on the sky and distances to stars based 
on a mass distribution model of the Galaxy. The stel- 
lar brightness is assigned based on a model luminosity 
function of stars considering distances to stars and ex- 
tinction. 

2. Once source stars are produced, we then produce lenses 
whose locations and masses are assigned based on mod- 
els of the mass distribution and mass function, respec- 
tively. 

3. To estimate the event rate in terms of events per year, 
it is required to model velocity distributions of source 
stars and lenses. The velocity distribution is also re- 
quired to estimate the distributions of the physical pa- 
rameters of lens systems and the observables of lensing 
events. 

The details of the simulation are described in the following 
subsections. 

2. 1 . Stellar Distribution 

We model the mass distribution in the solar neighborhood 
as a double-exponential disk of the form 



p(R,z) oc exp 



R-Ro^\z[ 



(1) 



where Rq is the galactocentric distance of the sun and /is and 
/i, are the radial and vertical scale heights, respectively. We 
adopt Ro = 8 kpc. For the radial scale height, we adopt a fixed 
value of hu = 4.0 kpc. For the vertical scale heights, on the 
other hand, we adopt varying values depending on the stellar 
brightness such that /i- = 117 pc for My < 2.0, 195 pc for 
Mv < 3.0, 260 pc for My < 4.0, 325 pc for My < 5.0, 390 pc 
for Mv < 6.0, 455 pc for My < 7.0, and 520 pc for My > 7.0. 

We model the luminosity function of stars by adopt- 
ing flia^_of__stos_Jnthe__solar neighborhood presented 
in iBinnev & Merrifieldl (119981) . This luminosity func- 
tion is constructed by com bining the data publishe d in 
Allen ('1973) for Mv < in iJahreiss & WielenI (Il983h and 
Kroupa, Tou t & Gilmord (|i990|) for My > 0. 

Extinction is modeled such that the stellar flux decreases 
exponentially with the increase of the dust column density, 
i.e. 

Av = -2.51og[exp(-/:I]d)], (2) 

where Ed is the dust column density and k is a proportional 
constant. The dust column density is computed on the basis 
of an expontial dust distribution model, i.e. 



Ds 
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Pd (X exp -- — , (3) 



where the integral is along the line of sight toward the source 
star and hd~ is the vertical scale height of the dust distribution. 
We adopt /ij = 150 pc. 

With these models, we produce source stars with their dis- 
tances from the sun and locations on the sky assigned based 
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Fig. 1. — Position of a star as seen from the sun. (Lb) are tlie galactic 
longitude and latitude of the star, d is the distance to the star, and So is the 
galactocentric distance of the sun. The vectors f , ip, and 6 represent the unit 
vectors in the spherical coordinates centered at the position of the sun. The 
projected velocitw of the stm' con'esBOfids to Vij =Viaii+f>g0. , 

on-'tne mass distribution model and their aDsomte magnitudes 
allocated based on the luminosity function. Then, the appar- 
ent magnitude is computed based on the distance to the star 
and the amount of extinction. We set the normalization of the 
stellar number density and extinction so that the star count 
result from our simulation matches the star count data as a 
function of galactic latitude presented in Allen ( 1998). 



2.2. Lens Distribution and Mass Function 

In our simulation, we assume that only stars are responsible 
for lensing events and no MACHOs or stellar remnants are 
taken into consideration. For the mass function of lens matter, 
we adopt the model of Gould (2000) that is constructed based 
on the observations of Zoc cali et al.. (.2000) . The model has a 
double power-law distribution of the form 



dN ^ / m 
dm \ ifir 
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-2.0 
-1.3 



for m > nic 
for m < Mr 



(4) 



where the turning-point mass is = 0.7 Mq. The distance 
to the individual lenses are assigned from the same mass dis- 
tribution model as that of source stars in equation ([Til. Since 
the lens is a star, its flux contributes to the observed flux. In 
this case, the flux from the lens works as blended flux to the 
observed flux. Then, the apparent m agnification of a n event 
appears to be lower than actual values (lNemiroflll997h . How- 
ever, considering that source stars are bright enough to be 
monitored with small telescopes while most lenses are faint 
stars, the lens c ontribution to the observed flux is in most 
cases negUgible (iHanll 19981) . 

2.3. Velocity Distribution 
We model the velocity distribution as gaussian, i.e. 

"(v-v)2" 



/(v) cx exp 



2^2 



(5) 



The adopted means and standard deviations of the distribu- 
tions of the individual velocity components in the galactocen- 
tric cylinderical coordinates are (v's , , v'j.) = (0, 220, 0) km s~' 
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and {crR,a^,a.) = (38,25,20) km s~', respectively, i.e. a flat- 
rotating disk with gaussian velocity dispersion. For lensing, 
only the projected velocity, vj^, as seen from the observer are 
related to lensing events. For lensing simulation, then, it is 
required to convert the velocity components into those in the 
spherical coordinates that are centered at the position of the 
sun. This conversion is done by the relation 

[ Vr \ / cos/? sin/? ^ 
10 
V -sin/j cos/j I 



cos a sin a \ / vr 
-sin a cos a 1 

1 / V V, 



(6) 

where b is the galactic latitude of the star and a is the angle 
between the two lines connecting the sun, the projected posi- 
tion of the star on the galactic plane, and the galactic center 
(see Figure[Tll. The angle a is related to the galactic longitude 
of the star (/), galactocentric distance of the sun, and the pro- 
jected distance to the star on the galactic plane from the sun 
(d) by 



a = sm 



Rq sin I 



= {rI + cF -2Rod cos I) 



1/2 



(7) 



where R^, is the projected distance to the star from the galactic 
center (see Figure[l]l. Then, the projected velocity of a star de- 
pends not only on the galactic coordinates {l,b) but also on the 
distance to the star d = dsecb, i.e. \j_ = {v,^,vg) = \±{l,b,d). 
With the projected velocities of the lens, Vl = \±{Di}, and 
source, Vl = Vx(Z)s), the relative lens-source transverse ve- 
locity is computed by 



Vt = Vl- 



, ( Ds-Dl 



(8) 



where vo represents the projected velocity of the observer and 
Dl and Ds are the distances to the lens and source, respec- 
tively. We note that the observer, lens, and source are all in 
rotation with the same rotation speed, and thus the contribu- 
tion of the rotation to the transverse speed is negligible. Then, 
the lens-source transverse motion is mostly caused by the dis- 
persion of the rotation-subtracted residual velocity. 

2.4. Event Rate 

Once the source, lens, and their relative speed are chosen, 
its contribution to the event rate is computed by 

r (xnil,b,Ds)n(l,b,DODlavt, (9) 

where n(l,b,d) represents the number density of stars at the 
position (l,b,d) and a is the cross-section of the lens-source 
encounter. Here the factor is included to account for the 
increase of the number of source stars with the increase of 
the distance. Under the definition of a lensing event as the 
approach of the source star within the Einstein ring of the lens 
(equivalently, source flux is magnified with a magnification 
A > 1.34), the lensing cross-section is set as cr = Ir^, where 
is the physical Einstein ring radius. The Einstein radius is 
related to the physical parameters of the lens system by 

^^ 1/2 



Dl(Z)s-£>l) 



Ds 



(10) 



where M is the mass of the lens. We assume that the survey is 
conducted all through the year' and events are detected with 
100% efficiency. 

' This requires a network of wide-field telescopes distributed over the 
Earth. 



TABLE 1 
Rates of Microlensing Events 



Vhm 


A'*,tot 


(r) 


r,ot (yr-1) 


12 


0.14 X 10^ 


0.16 X lO"** 


0.18 


14 


0.65 X 10^ 


0.32 X lO"** 


0.61 


16 


2.4 X 10^ 


0.60 X 10->* 


3.86 


18 


8.0 X 10' 


1.0 X 10-** 


21.20 



Note. — Rates of near-field microlensing 
events, Ftot. expected from all-sky surveys with 
various magnitude limits, Vjj,^. Also listed are the 
total number of stars that can be monitored from 
the surveys, W+.tot, and the average optical depth 
to lensing, (r). 

3. RESULTS 

In Table [T] we list the total event rate, Ftot, expected from 
surveys with various magnitude limits. Also listed in the table 
are the total number of stars that can be monitored from the 
surveys, A^*,tot, and the average optical depth to lensing, (t). 
The optical depth is determined from the mass distribution 
model by 



T{l,b)-- 



dDsn(l,b,Ds;V < Viim)DlJ^'dDLn(l,b,DOrl 



J^^dDsn(l,b,Ds;V <VundD\ 



(11) 



where the factor dDsn(l,b,Ds',V < Viim)D^ corresponds 
to the number of stars brighter than the magnitude limit 
toward the observational field, A^*, and the other factor 
Jq^^ dDi^n(l,b,Di)r^ is the accumulation of the area occupied 
by the Einstein rings of the individual lenses along the line 
of sight toward the source star We note that although there 
exists limitation in the brightness of source stars {V < Viim), 
there is no restriction in the lens brightness. The presented 
value in Table [T| are the mean optical depth averaged over 

180° 90° 

the whole sky, i.e. (t) = (4Tr)~^ J_^^Q„d I J_gQ^ T(l,b) sin b db, 
and the total number of stars over the entire sky, i.e. A^*,tot = 

180° 90° 

S-\iQidl J_gooNi,(l,b)smb db. We find that our estimation 
of the optical depth matches well with the estimation of 
[Nemiroff (1998). The multiplication A^^.tot x (t) represents 
the number of stars undergoing lensin g magnification at a 
given moment, which was presented in Nemiroff) (Il998l) . We 
also find a good match between his and our estimations. The 
estimated event rate ranges from Ftot ~ 0.2 per year for a sur- 
vey with a magnitude limit of Vnm = 12 to Ftot ^ 20 per year 
for a s urvey with Viitn =18, confirming the result of lNemirofB 
(Il998h that the lensing probability rapidly increases with the 
increase of the magnitude limits. Two factors contribute to the 
increase of the event rate with the increase of the magnitude 
limit. The first factor is the increase of the number of source 
stars and the other factor is the increase of the line of sight 
toward source stars and thus increase of the optical depth. We 
find that the former factor is more important. 

In Figure 111 we present the distributions of lens parameters 
and the observables of events. From the top, the individual 
panels represent the distributions of the source star brightness, 
distances to source stars and lenses, lens-source transverse 
speed, event time scale, and lens-source relative proper mo- 
tion, respectively. The event time scale is determined by the 
Einstein time scale, which is required for the source to tran- 
sit the Einstein radius of the lens, i.e. = te/v. The proper 
motion corresponds to /i = 9-E/t-E, where 0-^ = r^/D]^ is the an- 
gular Einstein radius. In Table |2l we also present the average 
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TABLE 2 

Average Values of Lens Parameters 



Mim (Os> (Dl> (v.) (fE> (M> 

(mag) (kpc) (kpc) (kms"') (days) (masyr"') 



12 0.39 0.19 55.3 18.7 71.5 

14 0.56 0.27 55.4 21.9 58.5 

16 0.78 0.36 55.6 25.0 49.1 

18 1.03 0.44 56.0 27.4 44.1 
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Fig. 2. — Distributions of lens parameters and observables of near-filed 
microlensing events expected from all-sky surveys with various magnitude 
limits, Viin,. From the top, the individual panels represent the distributions of 
the source star brightness (V), distances to the source star (Z)s) and lens (Dl), 
transverse speeds (vt), event time scale (?e), and lens-source relative proper 
motions (fi), respectively. The mean values of the parameters and observables 
are listed in Table|2] 

values of the lens parameters. From Figure |2] and Table [T] and 
|2]we find the following tendencies. 

1. The distributions of Ds and Dl vary considerably de- 
pending on the magnitude limit. The trend, as expected, 
is that the distances to the source and lens become 
larger as fainter stars are monitored. 

2. On the other hand, the dependency of the distribution 
of on the magnitude limit is not very strong and the 
dependency of the velocity distribution is nearly negli- 
gible. 

3. With the increase of the magnitude limit, the distance 



Note. — Average values of the lens parameters and 
observables of near-field microlensing events expected from 
wide-field surveys with various magnitude limits, Vijn,. 

to the lens increases while the transverse speed remains 
nearly the same. As a result, the average proper motion 
/i = 6'E/fE = Vt/Z^L of events decreases as the magnitude 
limit increases. However, we note that even for sur- 
veys with V\im ^18, the mean proper motion of events 
is (/i) > 40 mas yr"', which is much larger than the 
typical value of 5 mas yr~' of galactic bulge events. 
Then, the lens and source of a significant fraction of 
near-field events can be resolved from follow-up obser- 
vations by using high-resolution instrument such as the 
Hubble Space Telescope conducted several years after 
the peak of the magnification. This not only enables 
the measurement of the proper motion but also helps to 
identify the lens. 

Figure[3]shows the variation of the event characteristics and 
probability of lensing depending on the position of the sky. 
From the top, the panels in each column represent the dis- 
tributions of the average distances to source stars and lenses, 
lens-source transverse speed, event time scale, number den- 
sity of stars in the field, optical depth, and the event rate per 
unit angular area of the sky. The panels in each row represent 
the distributions for different magnitude limits. The tenden- 
cies found from the maps are as follows. 

1 . For a survey with a given magnitude limit, the average 
distances to source stars and lenses become smaller as 
the latitude of the field increases. This is because the 
number density of stars along the line of sight decreases 
more rapidly as the field is located further away from 
the disk plane. 

2. The transverse speed increases with the increase of the 
latitude. This tendency is due to the elongation of 
the velocity ellipsoid of the rotation-subtracted resid- 
ual stellar motion. The major axes of the velocity el- 
lipsoid are {(Tii,a^,a-) = (38,25,20) km s~' and thus 
the ellipsoid is elongated more along the galactic plane 
than along the pole direction. Then, the dispersion of 
the projected velocity distribution for stars located to- 
ward the direction of the galactic plane is (T±{b = 0°) = 
{a^ + al)^!^ 32 km s"', while the velocity dispersion 
for stars located toward the galactic pole direction is 
cfAb = 9Q°) = {crl + cfiyl^ - 45 km s"'. This also ex- 
plains the tendency of the longer time scales for events 
occurred on stars located at a lower latitude. 

3. With the increase of the stellar number density com- 
bined with the increase of the optical depth, the event 
rate increases rapidly toward the galactic plane. We find 
that > 50% of events are expected to be found in the 
fieled withZ7<20°. 
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Fig. 3. — Variation of the characteristics and probability of near-field microlensing events depending on the position of the sky. From the top, the panels in 
each column represent the distributions of the average distances to source star (Ds) and lens (Z)l), lens-source transverse speed (vt), event time scale (fg), average 
number density of stars in the field (W*), optical depth (r), and the event rate per unit angular area (F). The panels in each row represent the distributions for 
different magnitude limits. The map is centered centered at the galactic center, i.e. (Lb) = (0° ,0°). 



4. CONCLUSION 

Instigated by the recent discovery of an event occurred on a 
nearby star, we estimated the rate of near-field lensing events 
expected from all-sky surveys with various magnitude limits. 
Under the assumption that all lenses are composed of stars, 
our estimation of the event rate ranges from Ftot '^0.2 per year 
for a survey with a magnitude limit of = 12 to Ftot ^ 20 
per year for a survey with Viim =18, confirming the previous 
result that lensing probability rapidly increases with the in- 
crease of the magnitude limit. The increase of the rate is due 
to two factors which are the extension of the line of the sight 
toward source stars and the increase of the number of source 
stars. We found that the latter factor is more important. We 
also investigated the distributions of the physical parameters 
of lens systems and the observables of these events. From 
this, we found that although the average distances to source 
stars and lenses vary considerably depending on the magni- 
tude limit, the dependencies of the average event time scale 
and lens-source transverse speed are weak and nearly negli- 



gible, respectively. We also found that the the average lens- 
source proper motion of events expected from a survey with 
V\im =18 would be {fi) > 40 mas yr"' and the value further 
increases as the magnitude limit becomes lower. This implies 
that the source and lens of a typical near-field lensing event 
can be resolved from high-resolution follow-up observations 
conducted several years after the peak of the lensing magni- 
fication. From the investigation of the variation of the event 
characteristics depending on the position of the sky, we found 
that the average distances to source stars and lenses become 
shorter, the lens-source transverse speed increases, and the 
time scale becomes shorter as the the galactic latitude of the 
field increases. Due to the concentration of events near the 
galactic plane, we found that > 50% of events would be de- 
tected in the field with b <20°. 
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